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Abstract 

A solution of the old problem raised by Tolman, Ehrenfest, Podolsky 
and Wheeler, concerning the lack of attraction of two light pencils "moving 
parallel" , is proposed, considering that the light can be source of nonlinear 
gravitational waves corresponding (in the would be quantum theory of 
gravity) to spin— 1 massless particles. 

PACS 04.20. Jb, 02.40.-k, 04.30.-w 

Introduction 

In this paper we will study the repulsive behavior of gravitational interaction, 
in particle physics, associated to specific properties of some exact solutions of 
Einstein Equations. The interest in repulsive gravity, or antigravity as it was 
usually called, goes back to the fifty's [2 El E]. The general point of view was 
that since gravitational interaction is mediated by a spin-2 particle, it can only 
be attractive and thus, to obtain a repulsive behavior, some other ingredient 
is required. The idea was then to explore the possibility of repulsive matter- 
antimatter gravity, but within the old quantum field theories there was no room 
for such a possibility. The main arguments, reviewed in [3], were of various kinds 
including violation of energy conservation and disagreement with experiments 
of the Eotvos type due to the effects of antigravity on the vacuum polariza- 
tion diagrams of atoms. More recently however, within the context of modern 
quantum field theories, it was proven that those arguments were no longer suffi- 
cient to exclude repulsive effects and the interest in antigravity increased again. 
For example in [4] it was shown that in supergravity and string theory, due to 
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dimensional reduction, the effective 4-dimensional theory of gravity may show 
repulsive aspects because of the appearance of spin-1 graviphotons. 

Our point of view is the following. In the usual treatment of gravitational 
waves only Fourier expandable solutions of d'Alembert equation are considered; 
then it is possible to choose a special gauge (TT-gauge) which kills the spin-0 
and spin-1 components. However there exist (see section 2 and 3) physically 
meaningful solutions [5] [6] [7] [8] of Einstein equations which are not Fourier 
expandable and nevertheless whose associated energy is finite. For some of these 
solutions the standard analysis shows that spin-1 components cannot be killed 
[HI [TU] ; this implies that repulsive aspects of gravity are possible within pure 
General Relativity, i.e. without involving spurious modifications. In previous 
works it was shown that light is among possible sources of such spin-1 waves 

mug ma num. 

Photon-photon scattering can occur through the creation and annihilation of 
virtual electron-positron pairs and may even lead to collective photon phenom- 
ena. Photons also interact gravitationally but the gravitational scattering of 
light by light has been much less studied. Purely general relativistic treatments 
of electromagnetic wave interactions have been made resulting in exact solutions 
[T31 [2] , but these calculations are different from pure scattering processes and 
do not address the interaction at single photon level. It is not clear to what 
extent, calculations of the gravitational cross-section using QFT methods are 
consistent with classical GR. First studies go back to Tolman, Ehrenfest and 
Podolsky [15] and, later, to Wheeler [16] who analysed the gravitational field 
of light beams and the corresponding geodesies in the linear approximation of 
Einstein equations. They discovered that null rays behave differently according 
to whether they propagate parallel or antiparallel to a steady, long, straight 
beam of light, but they didn't provide a physical explanation of this fact. Later, 
Barker, Bathia and Gupta [17] , following a previous analysis of Barker, Gupta 
and Haracz |18j . analyzed in QED the photon-photon interaction through the 
creation and annihilation of a virtual graviton in the center-mass system; they 
found the interaction have eight times the "Newtonian" value plus a polariza- 
tion dependent repulsive contact interaction and also obtained the gravitational 
cross sections for various photon polarization states. Results of Tolman, Ehren- 
fest, Podolsky, Wheeler were clarified in part by Faraoni and Dumse [19], in the 
setting of classical pure General Relativity, using an approach based on a gen- 
eralization to null rays of the gravitoelectromagnetic Lorentz force of linearized 
gravity. They also extended the analysis to the realm of exact pp-wave solu- 
tions of the Einstein equations. After Barker, Bathia and Gupta, photon-photon 
scattering due to self-induced gravitational perturbations on a Minkowski back- 
ground has been analyzed by Brodin, Eriksson and Marklund [20 solving the 
Einstein-Maxwell system perturbatively to third order in the field amplitudes 
and confirming the dependence of differential gravitational cross section on the 
photon polarizations. 

Since the problem of the gravitational interaction of two photons is still 
unsolved, it appears necessary to take into full account the nonlinearity of Ein- 
stein's equations, just as in the case of gravitational waves generated by strong 
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sources [22 [22] • Section 1 contains a short account of Einstein field equations 
in the linear approximation, its gravitoelcctromagnetic formulation and corre- 
sponding application to the Tolman, Erhenfest, Podolsky problem. Section 2 
summarizes some geometric and physical properties (energy-momentum ten- 
sor and spin) of a family of exact solutions of Einstein equations, representing 
gravitational waves generated by a light beam or, more generally by massless 
particles. Section 3 is devoted to the motion of a massless spin-1 particle in the 
strong field regime. 

1 Weak Gravitational Fields 

1.1 The harmonic gauge 

A gravitational field g = g^ (x) dx^dx 1 ' is said to be locally weak if there exists 
a (harmonic) coordinates system and a region M' C M of space-time M in 
which the following conditions hold: 

flV = V* + V" << !> I << 1- (!) 

As it is known, in the weak field approximations in a harmonic coordinates 
system the Einstein field equations read 

□V = °- ( 2 ) 

The choice of the harmonic gauge plays a key role in deriving equation ([5]); no 
other special assumption either on the form or on the analytic properties of the 
perturbation h has been done. For globally square integrable solutions of the 
wave-equation (0) (that is, solutions which are square integrable on the whole 
of M), with a suitable gauge transformation preserving the harmonicity of the 
coordinate system and the "weak character" of the field, one can always kill the 
"spin-0" and "spin— 1" components of the gravitational waves. However, in the 
following we will meet some interesting solutions which do not belong to this 
class. 

1.2 Gravitoelectromagnetism 

A slightly different point of view, which is useful in clarifying the nature of spin 
of gravitational waves is provided by the gravitoelectromagnetism, henceforth 
GEM (see, for example, [25]). In this scheme one tries to exploit as much as 
possible the similarities between the Maxwell and the linearized Einstein equa- 
tions. To make this analogy evident it is enough to write a weak gravitational 
field fulfilling conditions (JTJ in the GEM fornfl 

ds 2 = c 2 (l + 2^-)dt 2 + -(A (9) • dx)dt - (1 - 2—=-)S ij dx t daf l , (3) 
d- c d- 

1 In this section the speed of light c will be explicitly written. 
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with 

"00 — ftOi — 5 — • 

C C 

Hereafter the spatial part of four- vectors will be denoted in bold and the stan- 
dard symbols of three-dimensional vector calculus will be adopted. In terms of 
and A^) the harmonic gauge condition reads 



1 

— + -V . A (9) 



and, once the gravitoelectric and gravitomagnetic fields are defined in terms of 
GEM potentials, as 

2c dt v ' 

one finds that the linearized Einstein equations resemble Maxwell equations. 
Consequently, being the dynamics fully encoded in Maxwell-like equations, the 
GEM formalism describes the physical effects of the vector part of the gravita- 
tional field. The situations which are usually described in this formalism are, 
typically, static: in fact, when this assumption is dropped, GEM gravitational 
waves are also possible. 

Then, the gravitoelectric and the gravitomagnetic components of the metric 
are given by 



where 

f$ = 0„4 ff) -M. fl) 

4 9) = -V/2= (-$ (9) ,A^). 

• The first order geodesic motion for a massive particle moving with velocity 

= (l,y), |v| << 1, in a light beam gravitational field characterized by 
gravitoelectric and gravitomagnetic fields, is determined (at first 
order in |v|) by the force: 

f(9) =_E (9) -2vAB (s) . 

• The first order geodesic motion for a massless particle moving with velocity 

= (1, v), |v| = 1, in the light beam gravitational field, parallel(anti) to 
z-axis (vj — iSjs) is slightly different 

f(s) = _ 2 +vAB^ 



The factor 2 in front of the gravitoelectric field can be understood as 
resulting from two contributions, one by the light beam, which is the 
source of gravity, and the other by the test photon [19]. It turns out 
p~5| [T9l [40] that for a massless particle moving parallel (antiparallel) to 
the light beam f^) = (f^ ^ 0). 
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2 Strong Gravitational Fields 



2.1 Geometrical properties 

In previous papers ( 26, 27, 28, S , [TQl [9] ) a family of exact solutions g of Einstein 
field equations, representing the gravitational wave generated by a beam of light, 
has been explicitly written 

g = 2f{dx 2 + dy 2 ) + n [(w (x, y) - 2q)dp 2 + 2dpdq] , (6) 

where p = A$ + B with A, B E K, $ (x, y) is a non constant harmonic func- 
tion, / = (V$) y/\p\/ p, and w (x,y) is solution of the Euler-Darboux-Poisson 
equation: 

Aw+ (d x In d x w + {d y In \p,\) d y w = p, 

A representing the Laplace operator in the (x, y) —plane and = pS^s^u3 the 
energy-momentum tensor. 

It is invariant for the non Abelian Lie agebra Q2 of Killing fields, generated 

by 

d d 

X = d- P > y = exp(p) a? 

with [X, Y] = Y, g (Y, Y) = and whose orthogonal distribution is integrable. 

In the particular case s = 1, / = 1/2 and fi = 1, the above family is 
locally diffeomorphic [9 to a subclass of Peres solutions [SJ [7] and, by using the 
transformation 

p = In |it| q = uv, 

can be written in the form 

g = dx 2 + dy 2 + 2dudv + -^du 2 , (7) 

with Aw(x,y) = p, and has the Lorentz invariant Kerr-Schild form: 
9iiv = f]jiv "i" Vkfj,kw> k^k^ = 0. 

2.2 Physical Properties 
2.2.1 Wave Character 

The wave character and the polarization of these gravitational fields can be 
analyzed in many ways. For example, we could use the Zel'manov criterion 
|29j to show that these are gravitational waves and the Landau-Lifshitz pseudo- 
tensor to find their propagation direction [51 |TD] . However, the algebraic Pirani 
criterion is easier to handle since it determines both the wave character of the 
solutions and the propagation direction at once. Moreover, it has been shown 
that, in the vacuum case, the two methods agree jlOj . To use this criterion, the 
Weyl scalars must be evaluated according to the Newmann-Penrose formulation 
|3T1 of Petrov classification [30] . 
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Thus, one has to choose a tetrad basis with two real null vector fields and 
two real spacelike (or two complex null) vector fields. Then, according to the 
Pirani criterion, if the metric belongs to type N of the Petrov classification, 
it is a gravitational wave propagating along one of the two real null vector 
fields. Since d u and d v are null real vector fields and d x and d y are spacelike 
real vector fields, the above set of coordinates is the right one to apply for the 
Pirani's criterion. 

Since the only nonvanishing components of the Riemann tensor, correspond- 
ing to the metric ©, are 

2 

Riuju = -jd i:j w(x, y), i, j = x,y 

these gravitational fields belong to Petrov type N [25]. Then, according to 
the Pirani's criterion, the metric ([7]) does indeed represent a gravitational wave 
propagating along the null vector field d u . 

It is well known that linearized gravitational waves can be characterized 
entirely in terms of the linearized and gauge invariant Weyl scalars. The non 
vanishing Weyl scalar of a typical spin— 2 gravitational wave is ^4. Metrics (J?} 
also have as non vanishing Weyl scalar ^4. 



2.2.2 Spin 

Besides being an exact solution of Einstein equations, the metric is, for 
w/u 2 << 1, also a solution of linearized Einstein equations, thus representing a 
perturbation of Minkowski metric 77 = dx 2 + dy 2 + 2dudv — dx 2 + dy 2 + dz 2 — dt 2 
(with u = (z — t)j\[2 v = (z + t)/\J~2) with the perturbation 

h := /loo = /133 = — ^03 = —^30 = To 

(z - tf 

generated by a light beam or by a photon wave packet moving along the z-axis. 

A transparent method to determine the spin of a gravitational wave is to 
look at its physical degrees of freedom, i. e. the components which contribute to 
the energy [32]. One should use the Landau-Lifshitz (pseudo)-tensor t£ which, 
in the asymptotically flat case, agrees with the Bondi flux at infinity [lOj . 

It is worth to remark that the canonical and the Landau-Lifchitz energy- 
momentum pseudo-tensors are true tensors for Lorentz transformations. Thus, 
any Lorentz transformation will preserve the form of these tensors; this allows 
to perform the analysis according to the Dirac procedure. A globally square 
integrable solution h^ v of the wave equation is a function of r = k^x* 1 with 
fc M /c M = 0. With the choice fc M = (1, 0, 0, —1), we get for the energy density tg 
and the energy momentum t$ the following result: 

16Trt = j( Ull -U22) 2 + u 2 2 , t° = t 3 

where = dh^ /dr. Thus, the physical components which contribute to the 
energy density are /in — /122 and h\i- Following the analysis of [32) . we see that 
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they are eigenvectors of the infinitesimal rotation generator 1Z, in the plane x—y, 
belonging to the eigenvalues ±2i. The components of which contribute to 
the energy thus correspond to spin— 2. 

In the case of the prototype of spin— 1 gravitational waves J7J), both Landau- 
Lifchitz energy-momentum pseudo-tensor and Bel- Robinson tensor [33J 133 135] 
single out the same wave components and we have: 

r o ~ c\{hox,x) 2 + C2(hoy :X ) 2 , t = t 

where c\ e ci constants, so that the physical components of the metric are h§ x 
and hoy. Following the previous analysis one can see that these two compo- 
nents are eigenvectors of ilZ belonging to the eigenvalues ±1. In other words, 
metrics ([7]), which are not pure gauge since the Riemann tensor is not vanish- 
ing, represent spin— 1 gravitational waves propagating along the z— axis at light 
velocity. 

Summarizing: globally square integrable spin—1 gravitational waves propa- 
gating on a flat background are always pure gauge. Spin— 1 gravitational waves 
which are not globally square integrable are not pure gauge. It is always possible 
to write metric ([7]) in an apparently transverse gauge [B|; however since these 
coordinates are no more harmonic this transformation is not compatible with 
the linearization procedure. 

What truly distinguishes spin—1 from spin— 2 gravitational waves is the fact 
that in the spin—1 case the Weyl scalar has a non trivial dependence on the 
transverse coordinates (x, y) due to the presence of the harmonic function. This 
could led to observable effects on length scales larger than the characteristic 
length scale where the harmonic function changes significantly. Indeed, the 
Weyl scalar enters in the geodesic deviation equation implying a non standard 
deformation of a ring of test particles breaking the invariance under of tt rotation 
around the propagation direction. Eventually, one can say that there should be 
distinguishable effects of spin—1 waves at suitably large length scales. 

It is also worth to stress that the results of [36l [3JJ [38] suggest that the 
sources of asymptotically flat PP_waves (which have been interpreted as spin—1 
gravitational waves (H [TO]) repel each other. Thus, in a field theoretical per- 
spective (see Appendix), pp-gravitons" must have spin—1 . 

3 Back to Tolman-Erhenfest-Podolsky problem 

In the previous strong field regime, geodesic motion can be still described in 
gravitoelectomagnetic formalism and we have 

Thus, "gravitational force" acting over a massless particle is given by 
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f(s) = -[w x (l - v z )i + w y (l - v z )j + (w x v x + w y Vy)k]/2u 2 . (8) 

Rather than geodesic orbits, the motion of spinning particles, should be 
described by Papapetrou equations 

n n etw i 

— { mv" + v (7 —) + -K fiU v*S^ = 0, 

where is the angular momentum tensor of the spinning particle and 
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S a = -e a ^ a vpS p<7 

defines the spin four-vector of the particle. These equations have been extended 
to the case of massless spinning particles by Mashhoon [231 125] ■ However, as- 
suming that the spin is directed along the z-axis S = (0,0, S z ), Papapetrou 
equations for photons coincide, in the gravitational field represented by Eq. ([7]) , 
with usual geodesic equations and no additional contributions must be added 
to the "gravitational force" given by Eq.®. 

The velocity v of a photon is determined by the null geodesies equations 

{h - 1) - 2hv z + (h + l)v 2 z = 

which has two solutions 

h — 1 w — u 2 

v z = 1, v z = 



h + 1 



If the photon propagates parallel to the light beam, v = (0,0, 1), then 

f (9) = 

and there is not attraction or repulsion. 

If the photon propagates antiparallel to the light beam, v = (h — 1) / (h + 1) . 
then 

f0?) =-Vw/2 (w + u 2 ) 
and the force turns out to be attractive. 



Conclusions 

Thus, the lack of attraction found by Tolman, Ehrenfest, Podolsky (later also 
analysed by Wheeler, Faraoni and Dumse) comes out also from the analysis of 
the geodesical motion of a massless spin-1 test particle in the strong gravitational 
field of the light, neglecting however the gravitational field generated by that 
particle. An exhaustive answer could derive only determining the gravitational 
field generated by two photons, each one generating spin— 1 gravitational waves. 
However, since helicity seems to play for photons the same role that charge 
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plays for charged particles, two photons with the same helicity should repel 
one another. This repulsion turns out to be very weak and cannot be certainly 
observed in laboratory but it could play a relevant role at cosmic scale and could 
give not trivial contributions to the dark energy. 

Therefore, together with gravitons (spin-2), one may postulate the existence 
of graviphotons (spin-1) and graviscalar (spin-0). Through coupling to fermions, 
they might give forces depending on the barion number. These fields might 
give [35] two (or more) Yukawa type terms of different signs, corresponding 
to repulsive graviphoton exchange and attractive graviscalar exchange (range 
>> 200m). 
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